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Abstract 

An ethylene trace gas technique was used to map out 
fluid transport and mixing within a circular-to-rectangular 
transition duct. Ethylene gas was injected at several 
points in a cross stream plane upstream of the transition 
duct. Ethylene concentration contours were determined at 
several cross stream measurement planes spaced axially 
within the duct. The flow involved a uniform inlet flow 
at a Mach number level of 0.5. Statistical analyses were 
used to quantitatively interpret the trace gas results. Also, 
trace gas data were considered along with aerodynamic 
and surface flow visualization results to ascertain transi- 
tion duct flow phenomena. Convection of wall boundary 
layer fluid by vortices produced regions of high total pres- 
sure loss in the duct. The physical extent of these high 
loss regions is governed by turbulent diffusion. 
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mean concentration defined by Eq. 

( 2 ) 

fluctuating concentration in Eq. (4) 
mean concentration coefficient 
defined by Eq. (3) 
duct inlet diameter 
turbulent diffusion coefficients 
(tensor notation) 

molecular diffusion coefficient in Eq. 

( 4 ) 

exponent in Eq. (1) 

radial distance from duct centerline 

radii in Eq. (1) 
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injection source strength in Eq. (7) 
mean velocity components 
mean velocity components (tensor 
notation) 

fluctuating velocity components 
fluctuating velocity components 
(tensor notation) in Eq. (4) 
Cartesian coordinates 
molar amount in Eq. (2) 
mean values defined by Eq. (12) 
variances defined by Eq. (13) 
covariance defined by Eq. (14) 


Introduction 

N ONAXISYMMETRIC exhaust nozzles are em- 
ployed on high performance aircraft to improve 
performance. For instance, the Pratt and Whitney two- 
dimensional/convergent-divergent nozzle described by 
Stevens et al } is a rectangular nozzle designed to al- 
low thrust vectoring, thrust reversal, jet area variation, 
and expansion area variation. Reduced thermal emis- 
sion is another benefit attributed to rectangular exhaust 
nozzles. Rectangular exhaust nozzles require a circular- 
to-rectangular transition duct to connect the engine and 
nozzle. 

Previous investigators have made detailed flow 
field surveys within circular-to-rectangular transition 
ducts using pressure probes, hotwire probes and laser 
anemometry. 2-7 These results indicate the dominant fea- 
ture within this flow field are pairs of vortices which ap- 
pear downstream of the transition region of the duct near 
the diverging sidewalls. Associated with these vortices 
are regions of substantial total pressure loss. 

For our investigation, a trace gas technique was used 
to determine fluid transport and mixing within a circular- 
to-rectangular transition duct. The trace gas technique 
involves injecting and tracking a nonreacting discernible 
foreign gas in the host flow field being investigated. In- 
formation about the host flow field is acquired by mea- 
suring the concentration of the trace gas in the host gas 
at many locations downstream of the point of injection. 
The trace gas technique is a valuable tool for several rea- 
sons; it is relatively inexpensive and easy to use, and it 
is a source of Lagrangian flow field data. 
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Our implementation of the trace gas technique in- 
volved injecting ethylene gas at seven different locations 
within a cross stream plane upstream of the transition 
duct and measuring the resulting flow field concentra- 
tion of ethylene at four cross stream measurement planes 
spaced axially within the duct. Twenty-eight different 
planer concentration distributions were measured. All 
measurements were made with uniform incoming flow 
and an inlet Mach number of 0.50. Some modifications 
were required to adapt the ethylene trace gas technique 
to the high speed flows involved in this study. Novel 
methods were employed to inject ethylene trace gas in 
the transition duct flow field with minimal disturbances 
and to sample the flow with virtually unlimited spatial 
resolution of the measurement locations in a cross stream 
plane within the transition duct. 

Statistical analyses were used to quantitatively in- 
terpret the trace gas results. In addition, the results of a 
detailed aerodynamic survey (e.g. flow velocity magni- 
tude and direction), free stream turbulence measurements 
and surface oil film visualization were considered with 
trace gas measurements to ascertain duct flow phenom- 
ena. Location, and to some extent the shape of the trace 
gas concentration contours is governed by convection of 
the trace gas by the mean flow. The amount of mixing 
or spreading of the contours is governed by turbulent dif- 
fusion of the trace gas. Convection of boundary layer 
fluid by vortices produced regions of lower total pressure 
levels (high loss). The physical extent of these high loss 
regions is governed by turbulent diffusion. 

Experimental Facilities 


Circular-to-Rectangular Transition Duct 

Fig. 1 shows the lower half of the circular-to- 
rcctangular transition duct and the Cartesian coordinate 
system used throughout this paper. The x-axis is coinci- 
dent with the duct centerline and the y- and 2 -axes are 
parallel to the major and minor axes of the duct exit. 
All coordinates have been normalized with the duct inlet 
diameter, which was 20.42 cm. In the yz-plane through 
each cross section, the surface of the duct satisfies Eq. (1) 
where the parameters r y , r z and n vary with the axial 
distance x and specify the exact cross-sectional shape. 



In the circular cross section region (n = 2), at 
x/D < 1.0, the duct maintains a constant circular cross 
section. The duct cross section is constant and nearly 
rectangular in the rectangular cross section region (n = 
10), at x/D > 2.5. The cross-sectional shape of the 
transition duct becomes more rectangular as the value of 
'* , '' ,, ?ases to a maximum value of 10, but it 


is never truly rectangular. The duct cross section changes 
from circular to nearly rectangular in the transition region, 
from 1.0 < x/D < 2.5. The aspect ratio of the duct, 
defined as the ratio of the major and minor axes lengths in 
the rectangular cross section region, is 3.0. The transition 
region Iength-to-diameter ratio is 1.5. 

The cross-sectional areas in the circular and rectan- 
gular cross section regions are equal. In the transition 
region the cross-sectional area increases to 1.15 times the 
area in the circular and rectangular cross section regions. 
The distribution of cross-sectional area in the transition 
region is as it would be if a duct with the same aspect 
ratio and transition region length-to-diameter ratio were 
constructed of circular and flat sections of metal only. 
This type of construction was believed to be representa- 
tive of fabrication methods used in the aviation industry. 

The duct is symmetric with respect to reflections 
across the horizontal xy- and vertical zz-planes. For flow 
without inlet swirl the flow field retained the duct’s two 
planes of symmetry. That is, any scaler flow field prop- 
erty, /, satisfied the condition f(x , y, z ) = f(x , -y, z) = 
f(x,y,-z) — /(x,-y, -z). Thus, measurements in 
only one quadrant of the duct were necessary. Measure- 
ments were made in two adjacent quadrants, however, in 
the upper half of the duct on either side of the zz-plane, 
to confirm that the measured flow was symmetric. 



Fig. 1 Circular-to-rectangular transition duct. 

Test Facility 

All transition duct experiments were conducted at 
the NASA Lewis Research Center using the Internal 
Fluid Mechanics Facility. This facility was designed to 
support tests of a variety of internal flow configurations. 
For the experiment reported here, air proceeded from 
the test cell through a large settling chamber containing 
a honeycomb and two screens, and thence through an 
axisymmetric contraction having an area ratio of 57:1. 
From the contraction, the flow passed through a straight 
pipe to allow a boundary layer to develop then into 
the transition duct. After passing through the transition 
duct the flow was exhausted into a discharge plenum 
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which was continuously evacuated by central compressor 
facilities. 


Trace Gas Measurement Equipment 


In this study, ethylene was used as the trace gas. 
Ethylene possesses several properties which make it de- 
sirable for trace gas studies: 

(1) The molecular weight of ethylene, 28.1, is nearly 
the same as air, 29.0. This eliminates the effect of 
buoyancy forces on the motion of the trace gas. 

(2) Under normal circumstances, ethylene and other hy- 
drocarbon gases are only found in very minute quan- 
tities in air. Thus, the problem of having to adjust 
the trace gas concentration measurements to com- 
pensate for the background level of trace gas in air 
is avoided. 

(3) The mean concentration of ethylene, or any hydro- 
carbon gas, can be measured very accurately with 
a flame ionization detector. This instrument is de- 
scribed later in this section. 

The equipment necessary to make ethylene trace gas 
measurements in an air flow can be divided into three 
major components; the ethylene trace gas injection sys- 
tem, the flow field sampling system and the ethylene con- 
centration measuring system. The relationship between 
these three major components and their subcomponents 
is shown in Fig. 2. Because there are no commercially 
available ethylene trace gas equipment items specifically 
intended for aerodynamic tests like those presently de- 
scribed, the equipment must be assembled from compo- 
nents designed for other applications. 

Ethylene Trace Gas Injection System 

The ethylene trace gas injection system was designed 
to introduce the ethylene gas into the flow field at a 
point location with a specific flow rate as unobtrusively 
as possible. The ethylene injection probe was designed 
to minimize the downstream flow disturbance from the 
probe while being sufficiently sturdy to withstand the 
aerodynamic forces of flows with Mach numbers as great 
as 0.5. Further, the probe had to be capable of injecting 
trace gas at an adequate rate to allow detection throughout 
the downstream flow field. 

The most common trace gas probe injection method 
previously used was an L-shaped probe, like a down- 
stream facing total pressure probe. Several studies that 
have employed the ethylene trace gas technique have in- 
cluded the results of experiments designed to measure 
and minimize disturbances from L-shaped type injection 
probes. 8-12 In these studies the diameter of the injection 
probes ranged from 0.56 mm to 3.2 mm. The results of 
these studies indicate that in regions of low turbulence 
intensity some minor probe disturbance is unavoidable 



(1) Ethylene supply 

(2) Mass flow controller 

(3) Injection probe 

(4) Sampling probe 

(5) Sampling pump 

(6) Total hydrocarbon 
analyzer 

(7) Hydrogen supply 

(8) Zero-gas air supply 

Flow field sampling system 


Ethylene trace gas injection system 
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Ethylene concentration measuring system 


Fig. 2 Trace gas measurement equipment schematic. 

and the smaller probes produced less trace gas spread- 
ing. However, in regions of higher turbulence intensity, 
the disturbance created by the probe is not of significant 
magnitude compared to the existing turbulence to mate- 
rially effect the mixing of the trace gas. 

All of the ethylene trace gas injection probes de- 
scribed in the cited literature were operated in compar- 
atively low speed flows and utilized a cantilever design. 
For the transition duct study presently considered a probe 
design change was necessary. For the injection probe to 
be as small as possible and yet withstand the aerody- 
namic loading imposed when it is extended to the cen- 
ter of the duct free stream at Mach numbers as great as 
0.50 required abandoning the cantilever design of previ- 
ously used probes. Instead, a catenary design was used as 
shown in Fig. 3. The injection tube was supported along 
its full length by a second tube which was held in tension 
across a rigid frame. The cross section of the rigid frame 
and the wall bounding the flow are shown in Fig. 3. The 
ethylene trace gas injection probe was constructed from 
stainless steel hypodermic needle tubing with an outside 
diameter of 0.56 mm. The catenary design allowed the 
ethylene trace gas injection probe to be as slender as the 
smallest probes reported in the cited literature, yet be op- 
erated at significantly higher flow field velocities than is 
capable with a cantilever configuration. Trace gas mea- 
surements in the transition duct indicate that the injection 
probe disturbance effects were minimal. 

The flow rate of ethylene trace gas was chosen so 
that the bulk velocity of ethylene at the injection probe 
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exit would approximately match the free stream veloc- 
ity. For the transition duct experiments, with the free 
stream Mach number nominally equal to 0.50, an ethy- 
lene volumetric flow rate of approximately 600 ml/min 

was required. The rate at which trace gas was injected 
was controlled with a Datametrics model 825 mass flow 
controller. 

Ethylene trace gas was injected from the seven dif- 
ferent locations shown in Fig. 1. All trace gas injection 
locations were within an axial plane situated one radius 
upstream of the beginning of the transition duct (x/D = 
— 0 . 5 ). One injection location was coincident with the 
centerline (y/D = 0 . 0 , z/D = 0 . 0 ). The remaining six 
injection locations were placed radially from the center- 
line at r/D = 0.249 and r/D = 0.498 along the z-axis 
(i.e. y/D = 0 . 0 , z/D = 0.249 and y/D = 0 . 0 , z/D = 
0 . 498 ) and along lines through the centerline 60 ° on either 
side of the z-axis (i.e. y/D = - 0 . 215 , z/D = 0 . 124 ; 
y/D = 0 . 215 , z/D = 0 . 124 ; y/D = - 0 . 430 , z/D = 
0.249 and y/D = 0 . 430 , z/D = 0 . 249 ). 


Flow Field Sampling System 

The purpose of the flow field sampling system was 
to provide the ethylene concentration measurement in- 
strument with a steady and continuous supply of flow 
field sample from a specified location. The components 
of the flow field sampling system included the sampling 
probe or tube, the probe positioning system, and the sam- 
ple pump and tubing necessary to maintain the flow of 
sample from the sampling probe to the ethylene concen- 
tration measurement instrument. 

The most common flow field sampling probe design 
used in the past employed features similar to a single-hole 
Pitot type pressure probe. 8, 9 * 13 The flow field sample 
was drawn through the single probe opening. Another 
probe design reported was similar to a three-hole "cobra” 
type pressure probe. 10 When using this kind of probe the 
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Fig. 3 Cantilever and catenary 
injection probe designs. 


two side holes were used to align the yaw axis of the 
probe with the flow field direction while the center hole 
was used to draw the sample. 

For the references cited above the measurements 
were performed at relatively low velocities, in all cases 
the Mach number was less than 0.10. One impediment to 
applying the ethylene trace gas technique to higher speed 
flows has been the perceived need for isokinetic sampling 
of the flow field (i.e. sampling at a velocity through the 
probe opening that matches the local host flow velocity). 
Isokinetic sampling is actually not necessary to accurately 
measure the local concentration of trace gas in a host 
gas flow. Injecting trace gas into the flow field creates 
an Eulerian concentration field. The local value of the 
concentration in the flow field does not depend on the rate 
at which the sample is acquired. However, if the sample 
velocity through the probe opening greatly exceeds the 
local host flow velocity the spatial resolution ability of 
the probe would deteriorate. 

Another sampling consideration is the importance of 
aligning the sampling probe precisely in the direction of 
the host flow. This concern is also largely unnecessary if 
the proper precaution is taken. Because of the measure- 
ment technique, the concentration of the sample will not 
depend on the orientation of the probe when the sample is 
acquired, provided the presence of the probe doesn’t ma- 
terially alter the flow field. The sampling probe is usually 
the same size as a pressure probe, where the disturbance 
caused by the probe is usually acceptable. These princi- 
pals of sampling rate and sampling probe alignment were 
confirmed during initial tests of the trace gas instruments 
and technique described in the Experimental Methods and 
Results section. 

The locations of the trace gas concentration mea- 
surement planes are shown in Fig. 1. The distance be- 
tween the probe access openings in the transition duct 
within each measurement plane was too large to satisfac- 
torily resolve the expected concentration gradients with a 
conventional flow field sampling probe. A novel method 
for sampling was therefore designed to allow finer spatial 
resolution in the measurement planes while still using the 
existing probe access openings. The design employed as 
the sampling probe a long straight flexible tube that ex- 
tended upstream to the sampling plane while firmly base 
mounted at a fixed location at the exit of the transition 
duct. The tube was constructed from stainless steel hypo- 
dermic needle tubing with outside and inside diameters of 
1.83 mm and 1.37 mm, respectively. The position of the 
tube opening was controlled by a positioning arm and the 
flow field sample was drawn through the tube opening. 

The trace gas sampling probe and positioning arm 
are shown in Fig. 4. The location where the probe is 
mounted at the duct exit lies downstream beyond the 
right border of this photograph. The position of the probe 
opening can be visualized as lying on the surface of an 
imaginary scribed sphere whose radius is equal to the 
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Fig. 4 Flow field sampling probe and positioning arm. 


sample passes through the detector. This makes it very 
important to maintain a constant flow rate of sample 
through the FID. Hardware within the total hydrocarbon 
analyzer maintained a constant specified flow rate into 
the FID regardless of the sample flow rate into the total 
hydrocarbon analyzer, provided, of course, the flow rate 
into the total hydrocarbon analyzer exceeded the desired 
flow rate to the FID. The fraction of sample fluid not 
directed to the FID was discharged to the atmosphere. 
Therefore, fluctuations in the flow rate of the trace gas 
sampling pump did not affect the flow rate into the FID. 
The suggested flow rate of sample fluid into the FID 
specified by the manufacturer, 25 ml/min, was used. 


probe length. Because the probe length was sufficiently 
large, the deviation of the surface of the scribed sphere 
from the measurement plane was negligible. Different 
length sampling tubes were used for each measurement 
plane. 

A sampling pump was required to supply the ethy- 
lene concentration measuring system with a continuous 
supply of flow field sample. This pump was necessary 
to overcome the pressure loss in the tubing between the 
sampling probe opening and the ethylene concentration 
measuring system and because the total pressure of the 
flow field was slightly less than atmospheric pressure. 
The sampling pump was located upstream of the ethy- 
lene concentration measuring system. A K.N.F. Neu- 
berger, Inc., Model N05ANI pump was used. This is 
an oil free pump designed specifically for contamination 
free pumping of gases. Placing the pump upstream rather 
than downstream of the ethylene concentration measuring 
equipment avoids operation at below atmospheric pres- 
sure which allows the possibility of atmospheric air leak- 
ing into the ethylene concentration measuring equipment 
and diluting the sample. 

Ethylene Concentration Measuring System 

A GOW-MAC Instrument Company Model 23-500 
Total Hydrocarbon Analyzer was used for this research 
project. This instrument utilized a flame ionization detec- 
tor (FID) for measuring hydrocarbon concentration. The 
FID consists of a small burner within which fuel and 
oxidizer are mixed and burned. Hydrogen and zero-gas 
air (air containing less than 0.5 ppm total hydrocarbons), 
were used as fuel and oxidizer, respectively, so there were 
no hydrocarbons present in this flame. When the sample 
was mixed with the flame, any hydrocarbons were burned 
and the carbon atoms became ionized. The carbon ions 
and electrons passed between two electrodes, decreasing 
the resistance between the electrodes and thus permitting 
an electric current to pass. This current was directly pro- 
portional to the amount of carbon ions present. 

The FID responds to both the amount of carbon 
atoms present in the sample, and the rate at which the 


Experimental Methods and Results 


Test Conditions 

An aerodynamic survey of the flow field one radius 
upstream of the transition duct inlet ( x/D = —0.5) was 
conducted to establish test conditions. The centerline 
Mach number was 0.50. The Reynolds number, based 
on the centerline velocity and the transition duct inlet 
diameter, was 2,086,000. The boundary layer thickness 
was approximately 4 percent of the inlet diameter of the 
transition duct and the boundary layer shape factor was 
1.42, confirming turbulent flow. 


Trace Gas Measurements 

The flame ionization detector was calibrated by 
preparing a mixture of ethylene and zero-gas air of known 
concentration and supplying this as a sample to the total 
hydrocarbon analyzer. This calibration was performed 
periodically throughout the period of trace gas testing. 
After the initial calibration, it was seldom necessary to 
make any further adjustments. 

The concentration measured in trace gas surveys was 
the mean ethylene mole fraction, given by Eq. (2). In 
Eq. (2), n clh y] cnc refers to the molar amount of ethylene 
present in a sample and n refers to the total molar amount 
of the quantity of matter in the sample. The results are 
presented as a mean ethylene concentration coefficient 
C* y where the local mean ethylene concentration has 
been normalized by the peak value of mean ethylene 
concentration measured for that particular distribution, as 
shown in Eq. (3). 


^ _ n ethylene 

°ethylene ~ 


( 2 ) 


c „ _ ^ethylene 
^'cthylcne.pcak 

Prior to conducting the transition duct experiments, 
preliminary trace gas measurements were made in a cali- 
bration jet to insure that the trace gas equipment was oper- 
ating properly. These measurements were also performed 
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to reveal information about the trace gas technique itself. 
A test was performed to determine the sensitivity of the 
ethylene concentration measurements to the ethylene in- 
jection rate. Ethylene was injected at the centerline of the 
calibration tunnel jet. Two separate cross stream sample 
surveys were conducted approximately 110 mm down- 
stream of the injection point. The rate of ethylene injec- 
tion varied between the two surveys while all other test 
conditions were unchanged. A comparison of the mea- 
sured concentration distributions presented as concentra- 
tion coefficients showed very little difference. Within 
some reasonable limitations, the concentration coefficient 
appears to be insensitive to the rate at which ethylene is 
injected into the host flow. 

To test the sensitivity of the ethylene concentration 
measurements to a mismatch between the sample velocity 
at the probe opening and the host flow velocity, two 
cross stream sample surveys were compared. For both 
surveys ethylene was injected at the centerline of the 
calibration jet and the flow was sampled in cross stream 
planes approximately 110 mm downstream of the point 
of injection. In one case the jet Mach number was 
approximately 0.10. and in the second case the jet Mach 
number was approximately 0.60. The rate at which 
the sample was drawn was held constant, insuring a 
significant velocity mismatch between these two cases. 
In both cases the sample velocity did not exceed the host 
flow velocity. The results showed little difference in the 
concentration coefficient, demonstrating that isokinetic 
sampling is not necessary. 

To test the sensitivity of the ethylene concentration 
measurements to the sampling probe alignment, an L- 
shaped type sampling probe was employed. This probe 
was constructed with the probe opening in line with the 
probe shaft axis so that the probe could be rotated, chang- 
ing the probe orientation, without changing the location 
of the probe opening. Thus, this probe could be inten- 
tionally misaligned from the host flow direction without 
changing the probe opening location. Ethylene was in- 
jected at the jet centerline and measurements were made 
downstream at a single location. The probe alignment 
was changed approximately 30° on either side of the host 
flow direction without any discernible change in the mea- 
sured sample concentration. This evidence suggests, that 
at least within this range of angles, concentration mea- 
surements are insensitive to misalignment of the probe 
with the host flow. 

The trace gas measurements in the transition duct 
were made following the preliminary tests described 
above. Twenty-eight separate planer concentration dis- 
tributions resulted from introducing ethylene at the seven 
injection locations and measuring the resulting concen- 
tration field at the four cross stream measurement planes. 
The results for each of the 28 separate planer concentra- 
tion distributions were obtained by recording the ethy- 
lene concentration at approximately 100 locations within 


the measurement plane. Ten separate ethylene concen- 
tration measurements were recorded at each location and 
the mean of these ten values was used for plotting and 
analysis. The 28 measured concentration distributions are 
shown in Fig. 5. The contour levels shown in Fig. 5 are 
at 20, 40, 60, and 80 percent of the peak value. 

The axial locations of the four ethylene concentra- 
tion measurement planes and the measurement planes for 
the aerodynamic measurements reported in Ref. 7 were 
not the same. This was because of the difference in hard- 
ware required for each method. Despite the different lo- 
cations, at three of the trace gas measurement planes it 
was possible to obtain aerodynamic data, either directly 
(trace gas plane 3 and aerodynamic plane 2 were at iden- 
tical locations, while trace gas plane 2 and aerodynamic 
plane 1 were very close) or by interpolation (trace gas 
plane 4 was between aerodynamic planes 3 and 4). Fig. 
5 also shows the transverse components of velocity super- 
imposed over the trace gas concentration contours. No 
aerodynamic measurements were available at trace gas 
plane 1. 


Discussion and Conclusions 


Eq. (4) governs the transport of trace gas in a 
turbulent incompressible flow. The left hand side of Eq. 
(4) represents the convection of the mean concentration 
of trace gas by the mean flow, the first term on the right 
hand side represents the molecular diffusion of the mean 
trace gas concentration and the second term on the right 
hand side represents the turbulent diffusion of the trace 
gas concentration. 


r 5£ = D d (dC\ d _fa) 

'Ox* dxi 


( 4 ) 


Applying mixing length theory for homogeneous turbu- 
lence, the turbulent diffusion tensor can be written in 
terms of a symmetric tensor. Further, the molecular dif- 
fusion term is much smaller than the turbulent diffusion 
term. Under these conditions, Eq. (4) can be expressed 
as Eq. (5). 


3C_ - JL(n 0 £l\ 

Ui 9 Xi " dxi V ij dxj) 


( 5 ) 


When the trace gas technique is used to infer infor- 
mation about the flow field there is an ambiguity problem. 
The ambiguity problem arises when experimental mea- 
surements, which represent a solution to Eq. (5) with 
a known boundary condition, are used to determine the 
coefficients of Eq. (5), £/* and D{j, which are unknown. 
Although the solutions to Eq. (5), with some restrictions 
on the coefficients and boundary conditions, are unique, 
Eq. (5) with different coefficients and identical bound- 
ary conditions can give rise to identical solutions. The 
ambiguity problem is exacerbated by the fact that the ex- 
perimental data are not known everywhere, but only at a 
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Fig. 5 Concentration contours and transverse velocity components. 


finite number of locations where measurements are made. 
This makes it difficult to interpret trace gas results with- 
out either making some assumptions which simplify Eq. 
(5), or supplementing the trace gas data with additional 
information, or both. The supplemental information may 
be experimental data or numerical results. 

The following analysis will consider one special case 
of the turbulent transport equation, when the molecular 
diffusion is negligible, the turbulent diffusion is homoge- 
neous but not necessarily isotropic, the turbulent diffusion 
in the axial direction is negligible, and the mean flow is 
uniform and in the x direction. Under these restrictions 
Eq. (5), an elliptic partial differential equation, reduces 
to Eq. (6), a parabolic partial differential equation like 
the heat equation, where the x dimension in Eq. (6) cor- 


responds to time in the heat equation. Without any loss 
in generality, it is assumed the y- and z - axes are coin- 
cident with the principal axes of the tensor D ijy so that 

Dyz ~ 0 . 


IT dC _ d 2 c d fdC' 

U " - Dyy dy* + 2 Dy *dy\dz, 


dx 


cfiC 


su 

_£/ 


(7) 

4nx^/Dy y D 22 

4x \ 

\ Dyy D l2 ) 

2 2xDj,j, 

ay= U 

= 

2 xD 22 

U 

(8) 
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Sampling 

Peak value location 

Mean value location 

Variance 

Covariance 

Area 

location 








estimate 

X 

D 

1 L 

D 

z 

D 

\L 1L 
D 

Mi 

D 

4 

4 

fry* 

D t 



(x 100) 

(x 100) 

(x 100) 

(x 100) 

(xlO 4 ) 

(xlO 4 ) 

(xlO 4 ) 

(xlO 4 ) 

0.97 

0.0 

0.1 

- 0.1 

0.4 

2.2 

3.2 

0.0 

8.3 

1.43 

0.0 

1.2 

- 0.2 

1.2 

2.8 

2.5 

0.4 

8.3 

1.99 

0.0 

0.7 

0.0 

0.4 

5.3 

1.8 

0.3 

9.6 

3.36 

0.1 

-0.1 

0.4 

0.2 

7.4 

3.0 

0.0 

14.8 


Table 1 Concentration distribution statistics for injection at y/D = 0.0, z/D = 0.0 


For the case of a point source of trace gas located 
at the origin, the solution to Eq. (6) is given by Eq. 
(7). The variances of the solution are given by Eq. (8), 
the covariance is zero because of the choice of axes. 
Although Eq. (6) appears simplified beyond the point of 
relevance, its solution provides some useful insight into 
interpreting trace gas results. For this particular solution: 

(1) The concentration contours will always be concen- 
tric ellipses. 

(2) The variances will always increase with axial dis- 
tance. 

(3) The area of the trace gas spreading can be estimated 
by Eq. (9) which is the area of the ellipse defined 
by Eq. (10). 


Area = - v\ z 


( 9 ) 


yV 2 


0-2 rj'i 

y z yz 


ZyZVyz 

0*0* — O 2 
u y u z u yz 


2 S 2 


0 * 0 * — o* 

y z yz 


= 1 (10) 


For each of the 28 measured concentration distribu- 
tions, the mean values, variances, and covariance of the 
concentration distributions were calculated. The defini- 
tions of these statistics are given in Eqs. (11) through 
(14). These values are presented in Tables 1 through 5. 
Also presented in these tables are the peak value location 


and an estimate of the spreading area based on Eq. (9). 
These statistics are useful to help quantify characteristics 
of the trace gas measurements. The variances measure the 
amount of trace gas spreading in the y and 2 directions 
while the covariance value indicates the orientation of the 
contours. All values in the tables are nondimcnsional. 


m = 


f 4>Cd,A 

JCdA ’ 


<t> = y>2>y 2 >y-z,2 2 


(ii) 


Hy = £[y], fi 2 = E[z ] (12) 


al = E 


(y — Vy) . v] - E (z- Hz) 


(13) 


<V = E[(y - Hy)(z - fi 2 )\ .(14) 

For trace gas injected at the transition duct center- 
line, y/D = 0.0, z/D = 0.0, the covariance, Table 
1, remained relatively small along the length of the duct. 
This indicates that the turbulence field is symmetric along 
the duct centerline. These trace gas contours became in- 
creasingly elongated in the y direction at downstream lo- 
cations. This can be seen in the magnitude of cr^/D 2 
relative to cr 2 /D 2 . Two possible explanations for this 


Sampling 

Peak value location 

Mean value location 

Variance 

Covariance 

Area 

location 








estimate 

X 

D 

JL 

D 

z 

D 

Mj l 
D 

Mi 

D 

4 

4 

fry t 



(x 100) 

(x 100) 

(x 100) 

(x 100) 

(xlO 4 ) 

(xlO 4 ) 

(xlO 4 ) 

(xlO 4 ) 

0.97 

- 0.1 

25.7 

- 0.2 

25.3 

2.1 

3.5 

0.2 

8.4 

1.43 

0.0 

23.9 

- 0.4 

24.3 

3.3 

2.5 

0.2 

9.0 

1.99 

0.0 

15.9 

0.0 

16.5 

5.6 

2.0 

-0.1 

10.4 

3.36 

1.5 

11.8 

0.5 

13.0 

9.9 

4.9 

-0.1 

21.8 


Table 2 Concentration distributions statistics for injection at y/D = 0.0, z/D — 0.249 
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Sampling 

Peak value location 

Mean value location 

Variance 

Covariance 

Area 

location 

X 

D 

1L 

D 

z 

D 

ESL 

D 

£jl 

D 


& 

fry* 

D 7 

estimate 


(xlOO) 

(xlOO) 

(xlOO) 

(xlOO) 

(xio 4 ) 

(xlO 4 ) 

(xio 4 ) 

(xlO 4 ) 


0.97 

-0.1 

49.4 

-0.1 

47.3 

6.2 

2.5 

-0.1 

12.3 

1.43 

0.0 

44.1 

-0.2 

42.4 

23.3 

2.3 

-0.3 

22.8 

1.99 

0.0 

28.4 

-0.4 

26.9 

40.3 

1.5 

0.1 

24.6 

3.36 

0.0 

25.0 

-0.3 

22.4 

50.7 

4.7 

0.0 

48.7 


Table 3 Concentration distributions statistics for injection at y/D = 0.0, z/D = 0.498 


are either the turbulence is anisotropic, that is v' 2 > w /2 , 
or the concentration distribution is being affected by the 
transverse components of velocity. Two experimental 
facts suggest that the elongation is produced by the trans- 
verse velocity components. First, the value of <j]/D 2 is 
decreasing from plane 1 to plane 3 and then increases 
from plane 3 to plane 4. It would be impossible for 
c 2 /D 2 to decrease by turbulent diffusion alone in a uni- 
form flow field. Second, velocity measurements show 
flow converging in the zz-plane and diverging in the xy- 
plane. The region of decreasing o\/D 2 corresponds to 
the region where the transverse velocity components are 
greatest. This can be seen in Fig. 5 at planes 2 and 
3. At plane 4 of Fig. 5, the transverse velocity compo- 
nents near the centerline are negligible and a 2 JD 2 has 
increased between planes 3 and 4. 

For trace gas injected at y/D = 0.0, z/D = 0.249 
(12 o’clock, halfway between the duct centerline and the 
wall) the statistical results, Table 2, arc similar to the re- 
sults for trace gas injected at y/D = 0.0, z/D — 0.0. 


However, at planes 3 and 4 the contour shapes for the 
two injection locations, Fig. 5, appear to be noticeably 
different. The contour shapes for trace gas injected at 
y/D = 0.0, z/D — 0.249 are no longer elliptic, rather 
they are bean shaped, having only one plane of symme- 
try. The simple model of homogeneous turbulence in 
a uniform mean flow allows only elliptic shaped con- 
tours, which would suggest that these contours result 
from turbulence inhomogeneity, or from nonuniformity 
in the mean flow. The mean flow certainly is nonuniform. 
However, there is nothing to indicate that the turbulence 
is not homogeneous. Thus, the bean shape of the con- 
tours probably results from transverse velocity compo- 
nents. Also, the amount of spreading is approximately the 
same as for trace gas injected at y/D — 0.0, z/D — 0.0, 
except at plane 4, which indicates roughly the same tur- 
bulence conditions in both regions of flow. 

For trace gas injected at y/D = 0.0, z/D = 0.498 
(12 o’clock, near the duct wall) the value of o\/D 2 in 
Table 3 indicates a very large elongation in the y direc- 


Sampling 

Peak value location 

Mean value location 

Variance 

Covariance 

Area 

location 

X 

D 

JL 

D 

z 

D 

ElL 

D 

El 

D 

4 

TP 

4 

fry* 

~d t 

estimate 


(xlOO) 

(x 100) 

(xlOO) 

(xlOO) 

(xlO 4 ) 

(xlO 4 ) 

(xlO 4 ) 

(xlO 4 ) 


-21.7 

12.1 

-21.4 

12.6 

3.9 

1.8 

-0.6 

8.2 

21.8 

12.5 

21.4 

12.4 

3.0 

3.0 

0.9 

8.9 

-25.8 

11.6 

-25.7 

12.1 

3.2 

2.6 

-0.6 

8.8 

25.6 

12.7 

25.6 

12.7 

3.0 

2.4 

0.8 

8.0 

-32.2 

9.7 

-32.0 

9.5 

5.3 

2.2 

-1.0 

10.3 

32.2 

9.4 

32.5 

9.4 

5.9 

1.8 

0.6 

10.2 

-40.8 

8.2 

-39.9 

8.2 

8.4 

3.3 

-1.4 

15.8 

38.3 

8.1 

38.2 

8.7 

8.2 

3.5 

1.7 

16.0 


Table 4 Concentration distributions statistics for injection at 
y/D = -0.215, z/D = 0.124 and y/D = 0.215, z/D = 0.124 
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Sampling 

Peak value location 

Mean value location 

Variance 

Covariance 

Area 

location 

X 

D 

JL 

D 

z 

D 

Ejl 

D 

D 

4 

w 

4 

O yt 

D t 

estimate 


(x 100) 

(xlOO) 

(xlOO) 

(XlOO) 

(xlO 4 ) 

(xlO 4 ) 

(xlO 4 ) 

(xlO 4 ) 



-42.1 

24.1 

-39.6 

23.4 

4.1 

5.5 

1.0 

14.7 

0.97 


42.2 

24.7 

40.3 

23.2 

3.9 

4.8 

-0.8 

13.4 


-46.6 

21.9 

-46.6 

23.0 

6.9 

7.1 

2.6 

20.5 

1.43 


49.9 

21.6 

46.3 

22.9 

6.9 

6.9 

-2.4 

20.3 


-62.4 

17.1 

-61.2 

15.9 

17.8 

25.9 

12.4 

55.0 

1.99 


62.5 

17.2 

62.9 

14.5 

17.3 

25.1 

-13.1 

50.4 


-72.7 

12.8 

-69.0 

9.7 

24.1 

55.1 

-5.4 

113.2 

3.36 


72.6 

14.9 

69.1 

9.9 

25.0 

53.8 

5.1 

114.2 



Table 5 Concentration distributions statistics for injection 

at 



y/D 

= -0.430, 

z/D = 

0.249 and y/D 

= 0.430, 

z/D = 

= 0.249 



tion. The results of Klebanoff 14 show that the magnitude 
of the transverse turbulent velocity fluctuations, and tur- 
bulent viscosity coefficient, increase greatly in a turbulent 
boundary layer. In a turbulent boundary layer the turbu- 
lence is very heterogeneous. The turbulence decreases 
rapidly away from the duct surface, in this case the z 
direction, while retaining its magnitude in directions par- 
allel to the wall, specifically the y direction. This should 
lead to preferential diffusion in the y direction. However, 
turbulent diffusion alone is not exclusively responsible 
for the observed contours because a l/D 2 again decreases 
from plane 1 to plane 3. The total amount of spreading 
is roughly twice that observed of trace gas injected at 
y/D = 0.0, z/D = 0.0 and y/D = 0.0, z/D = 0.249. 
Surface oil film visualization results reported in Ref. 7 
show large transverse velocities in this region. It appears 
that both turbulent diffusion and transverse velocities are 
responsible for the large transport of trace gas in the y 
direction displayed by the contours. 

For trace gas injected at y/D = -0.215, z/D — 
0.124 and y/D = 0.215, z/D = 0.124 (about 10 and 2 
o’clock, between the duct centerline and wall) the amount 
of spreading at each of the four planes is comparable to 
the spreading of trace gas injected at y/D = 0.0, z/D — 
0.0 and y/D = 0.0, z/D = 0.249. This indicates that 
outside the boundary layer the turbulence is nominally 
homogeneous. Each contour becomes more elongated 
as it travels downstream. Unlike the results for trace 
gas injected at y/D — 0.0, z/D = 0.0 and y/D = 
0.0, z/D = 0.249, the contours change their orientation. 
This can be observed in the covariance values. The major 
axes of the contours are no longer parallel to the y-axis, 
but are inclined to the y-axis. This probably results from 
the transverse velocities. 

Trace gas injected at y/D = -0.430, z/D = 0.249 


and y/D = 0.430, z/D — 0.249 (about 10 o’clock and 
2 o’clock, near the duct wall) produced the most interest- 
ing results. At planes 1 and 2 the statistics, Table 5, are 
nearly identical to the statistics for trace gas injected at 
y/D = 0.0, z/D = 0.498, except for the orientation. At 
plane 2, however, the contours in Fig. 5 appear visibly 
different, assuming a cusped teardrop shape. At plane 3 
the statistics differ significantly from the statistics of trace 
gas injected at y/D = 0.0, z/D — 0.498. The spreading 
was approximately twice as great. This resulted from the 
transfer by axial vortices of boundary layer fluid, with its 
greater turbulence intensity and diffusion, to a large area 
near the side walls. At plane 4 the much greater spreading 
by vortices is obvious. The unusual contour shape results 
from pairs of counter-rotating vortices near the transition 
duct side walls. These vortices were identified and dis- 
cussed in detail in Ref. 7. The total pressure coefficient 
contours of Fig. 6 clearly indicate large regions of lower 
total pressure near the duct side walls. The data of Figs. 
5 and 6, along with the results of the statistical analyses, 
confirm that the convection of boundary layer fluid by the 
sidewall vortices and turbulent diffusion are responsible 
for these large regions of lower total pressure. 

Summary 

The ethylene trace gas technique provided an accu- 
rate Lagrangian flow field representation for high velocity 
flow within the transition duct. These data are useful for 
assessing the accuracy of CFD codes used to predict the 
transport and mixing of a nonreacting gas. The methods 
developed to obtain the ethylene trace gas data can be 
used in other experimental configurations with very little 
modification. In fact, this technique is presently being 
used by Davis et al. is to measure the mixing of two su- 
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X 



Plane 4 
(x/D = 3.93) 


Plane 3 

(x/D = 2.55) 


Plane 2 
(x/D = 1.99) 


Plane 1 
(x/D = 1.49) 


Fig. 6 Total pressure coefficient contours. 


pcrsonic streams in a study of mixing nozzles for scramjct 
applications. The statistical methods used to help inter- 
pret the trace gas data provided a reliable way to quantify 
and compare contour characteristics such as shape, ori- 
entation and the amount of contour spreading. Trace gas 
data, along with aerodynamic and surface flow visualiza- 
tion data, confirmed the strong convection of wall bound- 
ary layer fluid by vortices to form significant regions of 
high total pressure loss in the duct. 
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